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The propagation of (massless) scalar, electromagnetic and gravitational waves on fixed 
Schwarzschild background spacetime is described by the general time-dependent Regge-Wheeler 
equation. We transform this wave equation to usual Schwarzschild, Eddington-Finkelstein, Painleve- 
Gullstrand and Kruskal-Szekeres coordinates. In the first three cases, but not in the last one, it is 
possible to separate a harmonic time-dependence. Then the resulting radial equations belong to the 
class of confluent Heun equations, i.e., we can identify one irregular and two regular singularities. 
Using the generalized Riemann scheme we collect properties of all the singular points and construct 
analytic (local) solutions in terms of the standard confluent Heun function HeunC, Frobenius and 
asymptotic Thome series. We study the Eddington-Finkelstein case in detail and obtain a solution 
that is regular at the black hole horizon. This solution satisfies causal boundary conditions, i.e., it 
describes purely ingoing radiation at r = 2M. To construct solutions on the entire open interval 
r € ]0,oo[, we give an analytic continuation of local solutions around the horizon. Black hole 
scattering and quasi-normal modes are briefly considered as possible applications and we use semi- 
analytically calculated graybody factors together with the Damour-Ruffini method to reconstruct 
the power spectrum of Hawking radiation emitted by the black hole. 


I. INTRODUCTION 

While many features of black hole spacetimes can be 
investigated in terms of the motion of (massive and mass¬ 
less) particles, wave propagation opens a perspective on 
additional phenomena such as interference effects, scat¬ 
tering of radiation at a black hole, its quasi-normal modes 
and black hole evaporation. Scattering and interference 
patterns have proven to be very important tools in other 
branches of physics and yield observational properties of 
different physical systems. Therefore, it is of great rel¬ 
evance to consider black-hole scattering as well and to 
describe the dynamics of fields in the presence of a black 
hole for different boundary and initial conditions. 

The central equation that describes the propagation of 
massless waves on Schwarzschild background spacetime 
is the general (time-dependent) Regge-Wheeler equa¬ 
tion which is a wave equation written in the usual 
Schwarzschild time and the so-called tortoise coordinate. 
Regge and Wheeler found this equation for spin s = 2 by 
investigating the stability of a Schwarzschild black hole 
under gravitational perturbations in a linearized pertur¬ 
bation theory [T]. To study electromagnetic wave prop¬ 
agation one can solve Maxwell’s equations on flat space- 
time but consider a medium with specified constitutive 
relations Hig. Another possibility to derive the Regge- 
Wheeler equation for s = 1 is to use the so-called Debye 
equation, because the source-free Maxwell equations re¬ 
duce to this single scalar equation on a Schwarzschild 
background n [5]. For scalar perturbations and their 
propagation the Klein-Gordon equation has to be em¬ 
ployed and the Regge-Wheeler equation for s = 0 fol- 
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lows from a straight-forward calculation [^. For each of 
the three different types of perturbations the respective 
Regge-Wheeler-type equation follows after an expansion 
of the fields into (tensorial, vectorial and scalar) spher¬ 
ical harmonics. Assuming additionally a harmonic time 
dependence, i.e., Fourier expanding the solution of this 
wave equation, leads to a radial equation that is known 
as the stationary Regge-Wheeler equation. This radial 
equation is a Schrodinger-type differential equation with 
a spin-dependent potential barrier. Due to the form of 
this potential no exact solution is known as long as the 
tortoise coordinate is kept as the radial variable. 

Changing the coordinate system and reconsidering the 
previous wave equations gives rise to new radial equations 
that admit analytic solutions in terms of series of special 
functions such as hypergeometric or Coulomb wave func¬ 
tions, see e.g. Ref. [7]. Another very insightful method 
is to use the Heun functions and the theoretical 

framework of singularity analysis to construct analytic 
(local) solutions to the radial equations. We will follow 
this strategy to a great extent in this article. 

We start with some necessary preparations in Sec. [l^ 
where we introduce the confluent Heun equation and 
important notations concerning singularities of differen¬ 
tial equations and local solutions around them. Fur¬ 
thermore, we consider on a Schwarzschild background 
a wave equation in the tortoise coordinate with an ar¬ 
bitrary potential term. This equation is transformed to 
usual Schwarzschild coordinates. By considering a gen¬ 
eral coordinate transformation that keeps the angles un¬ 
changed we proceed to Eddington-Finkelstein, Painleve- 
Cullstrand and Kruskal-Szekeres coordinates. In the last 
part of this section we sketch how to derive the respec¬ 
tive wave equations for the three different spin values. We 
combine the results to obtain the general Regge-Wheeler 
equation that has the form of the previously discussed 
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wave equation in the tortoise coordinate and can, thus, 
be transformed to the respective wave equations in the 
other coordinate systems. 

In Sec. |III| we study the radial equation in the tortoise 
coordinate, which is the stationary Regge-Wheeler equa¬ 
tion. A lot about the phenomenology can be learned from 
this simplest form of the radial equation. We describe 
asymptotic solutions that fulfill the causal boundary con¬ 
dition at the horizon, the so-called IN-mode for scattering 
and the QNM-asymptotics. Here and in the following, 
QNM stands for quasi-normal modes. Since there is no 
known exact analytic solution of the stationary Regge- 
Wheeler equation beyond the asymptotic behavior, ap¬ 
proximations or numerical methods have to be used to 
obtain the full solutions of this equation. Afterwards, we 
derive the radial equation in the usual Schwarzschild co¬ 
ordinate, which was already considered e.g. by Leaver 
[3 Ho] and Fiziev [H [9]. We briefly summarize their 
results and give properties as well as local solutions of 
this radial equation. All these solutions have only finite 
and non-sufficient domains of convergence and there is 
no regular solution at the black hole horizon that allows 
to penetrate this radius. 

In Sec. |IV| we derive and discuss the radial equation 
in Eddington-Finkelstein coordinates in full detail, and 
we examine its properties and local solutions in terms 
of Frobenius and Thome series and using the confluent 
Heun function HeunC. As far as we know, there is no pre¬ 
vious work concerned with the analysis of the Eddington- 
Finkelstein radial equation to find analytic solutions of 
the perturbation equations on Schwarzschild background, 
apart from a short remark on their asymptotic behavior 
at the horizon in Ref. m- To overcome the fact that 
the region where the local solutions converge is bounded, 
we show how an analytic continuation of the local Frobe¬ 
nius solutions onto the entire interval r s]0,oo[ can be 
performed. In this way we construct an exact solution on 
this interval that is regular at the horizon and satisfies 
the causal boundary condition. 

In Sec. jVjwe examine some applications of our analyti¬ 
cally continued solutions. We consider black hole scatter¬ 
ing and briefly sketch how to calculate QNM-frequencies 
for the different perturbations. At the end of this section 
we use the method of Damour and Ruffini m and our 
semi-analytically calculated graybody factors to reobtain 
the power spectrum of the Hawking radiation emitted 
from a Schwarzschild black hole. 


II. PREPARATION 


with rational coefficient functions Pi{z), where z shall be 
an element of the Riemann sphere, i.e., of the extended 
complex plane including z = oo. Following |12| . a point 
z = Zo 

an ordinary point if Vf Pi{z) is analytic a.t z = zq 
a singularity otherwise. 

A singularity at the point z = Zj 

regular if Vj {z — Zj)'^~'^ pi{z) is analytic at z = Zj 
irregular otherwise. 

If a differential equation possesses only regular singular¬ 
ities it is referred to as being Fuchsian, otherwise it is 
called a confluent equation. The singly confluent Heun 
equation (CHE) belongs to the remarkably large class of 
Heun’s differential equations. All members of this class 
can be derived from a single second order Fuchsian differ¬ 
ential equation of the form Q with four singularities, the 
general Heun equation (GHE), see e.g. [T^HTB] . There are 
four different confluent cases which are obtained through 
(consecutive) confluence processes of regular singulari¬ 
ties. In the case of the CHE, two regular singularities 
of the GHE, one of them located at infinity, merge and 
form an irregular singularity. One standard form of this 
CHE is given by Q, where the coefficient functions are 

Pi(^) = - H-(2a) 

z z — I 

P.(.) = - + 1 ) . (2b) 

and regular singularities are located at z = 0,1. By con¬ 
sidering the transformation C = 1/z we find an irregular 
singularity at C = 0: he., z = oo. The three singulari¬ 
ties can be further characterized by their singular rank 
(s-rank). The s-rank of a regular singularity is by defini¬ 
tion always equal to one, so we only have to calculate the 
rank of the irregular singularity. For the equation in the 
form ([^ this s-rank is found to be two, where we follow 
the convention in Ref. m- The generalized Riemann 
scheme (GRS), in the form proposed by Slavyanov and 
Lay [m, is a very useful tool in the analysis of singular 
differential equations. It contains 


A. The confluent Heun equation and its properties 


- s-ranks and locations of all singularities (1®* and 
2'^'^ row). 


We consider an ordinary, second order, linear and ho¬ 
mogeneous differential equation of the form 


- indicial (characteristic) exponents of the regular 
singularities (3“''^ and 4*^ row) and 


dz^ 


+ Pi(z) 


_d 


+ Po{z) 


y{z) = 0 


( 1 ) 


- characteristic exponents of second kind for the ir¬ 
regular singularity (3'^'^ row and following). 
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The corresponding GRS for the confluent Heun equation 
in the form ([^ is 

/II 2 \ 

0 1 oo ; z 

0 0 a ;q , 

1 — 7I — (57 + (5 — a 

0 

V /? / 

and yields all necessary information about properties of 
local solutions at the singularities of the differential equa¬ 
tion: Around both regular singularities local Frobenius- 
type solutions m can be constructed using the entries of 
the first and second column of the GRS and we can con¬ 
struct Thome-type solutions (asymptotic series) at the 
irregular singularity with the help of the third column. 


Frobenius solutions 

The two local Frobenius solutions around the regular 
singularity at z = 0 are 


00 

y^{z-,0) ^'^akz’^, (4a) 

00 

(4b) 

fc =0 

and those around the other regular singularity at z = 1 
are given by 

00 

y^{z'A)='^Ck{z-l)'', (5a) 

k^O 

00 

/^(z;l) =^4(z-1)'=+^"^. (5b) 

k=0 

The series expansion coefficients ak,bk, Ck and dk can be 
obtained by inserting the respective solutions from above 
into the GFIE § and solving the resulting three term re¬ 
currence relations with chosen initial conditions. Here 
and in what follows, we use the following notation: The 
superscript is either I or II to distinguish the two in¬ 
dependent local solutions, while the corresponding sin¬ 
gularity, at which the solution is constructed, is given 
in the argument after the semicolon. The local Frobe¬ 
nius solutions 2 /^’^^(z;-) are convergent within a circle 
in the complex plane, centered at the respective regular 
singularity with a radius that is the distance to the next 
neighboring singular point. Hence, in each case the re¬ 
gion of convergence is bounded by the unit circle centered 
at the considered singularity. We emphasize that at each 
regular singularity z = Zj S { 0 , 1 } the solutions have the 
asymptotic behavior 

~ (z — Zj)° = const, or ~ (z — Zj)^ for z ^ Zj , 


where the (complex) exponent A is either 1 — 7 or 1 — 5, 
depending on which of the two singularities is consid¬ 
ered. Another very useful canonical form of the CHE is 
obtained by choosing the coefficient functions in 0 to 
be 


pi(z) = 0 - 1 - 


b+1 

z 


c-l- 1 
z- 1 ’ 


Po{z) 



( 6 ) 


and, additionally, defining the two parameters d and e 
by 

fj, = -{a — b — c—2e + ab—bc) , (7a) 

iz = - {a + b + c + 2d + 2e + ac + be) . (7b) 

This form will be called the Maple-form of the CHE in 
the following, since it is implemented in the computer 
algebra system Maple [TB] . The GRS for the differential 
equation in the form ® then changes to 

/II 2 X 

01 00 ;z 

0 0 ^ 

-b -c b + c + 2-i^ 

0 

V -a / 


and yields the relations between the parameter sets 
{a,/?, 7 , i5, g} and {o, 6 , c, z/}. At a first glimpse, this 
looks like an unnecessary modification but it will prove 
very useful in the following. We define the standard con¬ 
fluent Heun function HeunC(a, b, c, d, e, z) in the domain 
|z| < 1 as the first local Erobenius solution at the regular 
singularity z = 0 [TB] 

00 

HeunC(a, 6 , c, d, e, z) := y^{z; 0) = ^ Ofc z^ (9) 

k=0 

with the additional normalization 


HeunC(a, 6, c, d, e, z)z=o =1; 


dHeunC(a, 6 , c, d, e, z) 


z=0 


5 -j- (c — a)(5 4- 1 ) 4- 2 e 

2 ( 6 + 1 ) 


Erom 0 we obtain the series coefficients ak as solutions 
of the recurrence relation 


Ofc = 0 Vfc < 0 (11a) 

Oo = 1 ( 11 b) 

/Ifc ^k tlk — 1 T ^k ttk—2 ; (11^) 

depending on the five parameters a, 6 , c, d and e, where 


Ak — 1 + 


6 

k ’ 


Bk — 1 — 


: — 6 — I 





( 12 a) 
, ( 12 b) 
( 12 c) 
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We can read off the behavior for large fc, since the limit 
A: —)■ oo yields 

lim Ak = lim Bk = 1, lim Ck = 0 

k—>-OC) fc—>-oo /c—>-oo 

^ Ofe — Uk-i —>■ 0 . 


Hence, we conclude that the standard solution 
HeunC(a, b, c, d, e, z) is convergent and, thus, well defined 
only for Iz| < 1 as it should be for a local Frobenius solu¬ 
tion; the radius of convergence is given by the distance to 
the next neighboring singularity. The major advantage 
is now that once the standard function HeunC is imple¬ 
mented together with the recurrence relation (12), the 


second Frobenius solution at z = 0 and even those two 
at the other regular singularity z = 1 can be expressed 
using this single function m- The second Frobenius so¬ 
lution at z = 0 can be represented as 


y \z;0) = z “HeunC(a, —6,c,d, e,z) 


(14) 


and the two Frobenius solutions at z = 1 are now given 
by 

y^{z; 1) = HeunC(—a, +c, b, —d, e -I- d, 1 — z), 

(15a) 

y^^{z; 1) = (z — 1)“'^ HeunC(—a, —c, b, —d, e -I- d, 1 — z). 

(15b) 

We will need these representations to construct solutions 
to the radial part of wave equations on Schwarzschild 
background spacetime in sections |III| and |IVl 


Thome solutions 

The Thome solutions (asymptotic series) at the irreg¬ 
ular singularity z = oo are constructed as m 

oo 

/(z;oo)= (16a) 

k^O 

oo 

y^\z] oo) = e^^ ^ au z-('=+t'+‘ 5-“) (16b) 

k=0 

and the coefficients pk and cr^ follow, again, from recur¬ 
sion relations that can be derived by inserting y^’^^{z; oo) 
into the differential equation. 


B. Wave equations in Schwarzschild spacetime 

For the entire article we choose natural units such 
that the speed of light, Newton’s Gravitational constant, 
Planck’s constant and Boltzmann’s constant are set to 
unity, i.e., c = G = h = ks = ^- Furthermore, all 
distances are measured in multiples of the Schwarzschild 
radius Vg = 2M. Using this convention, the horizon is 
located at r = 1 and the light sphere at r = 3/2, where 


r is the usual Schwarzschild radial coordinate (i.e., the 
areal radius). The Schwarzschild metric in these units is 
then 

guydx^dx’' = —— —^dt^ -|- -— - —-dr^ + r^dil^ , 
r [r -1) 

with dn^ being the standard metric on the unit two- 
sphere. 


Tortoise coordinates 


Since it will become very useful later, we start with an 
analysis of the wave equation 


drl 


dt"^ 


V{r) 


4'(t,r*) = 0 


(17) 


that involves the Schwarzschild time coordinate t and the 
tortoise coordinate 


r* = r -I- log(r - 1), (18a) 

dr* = dr = -—,ydr . (18b) 

dr (r — 1) 

Here V (r) is some given spherically symmetric potential. 
For U(r) = 0 we obtain the free wave equation, which is 
solved by 

4 -(t, r*) = A{w) -h H(w) , (19) 

i.e., by a superposition of in- and outgoing waves with 
arbitrary complex amplitudes A{oj) and B{uj). In the 
following we transform the wave equation 0 into the 
corresponding wave equations in usual Schwarzschild, 
Eddington-Finkelstein (EF), Painleve-Gullstrand (PG) 
and Kruskal-Szekeres (KS) coordinates. 


Schwarzschild coordinates 


To obtain the wave equation in usual Schwarzschild 
coordinates we perform the simple transformation r* i—> r 
that does not affect any of the other coordinates. Using 
(18) and the chain rule we rewrite 


d dr 9 (r — 1) (9 

dvi, dr* dr r dr 

(r - 1) 02 Id' 

r dr^ r^ dr 
and the wave equation 0 becomes now 


(r- 1) 



(20a) 

(20b) 


r — 1\ ^ 9^ {r — 1) d 
r ) dr'^ dr 


dt^ 


V{r) 


4'(t,r) 


0 . ( 21 ) 
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This equation contains a “new” first-derivative term and 
we notice that the values r = 0 and r = 1, i.e., the origin 
and the horizon in our coordinates, describe singularities 
of this differential equation. Starting from (21) we use a 
further coordinate transformation (t, r) i—>■ (u/v), where 


v = f{r,t), u = g{r,t), 


( 22 ) 


Eddington-Fmkelstein coordinates 

For (ingoing) Eddington-Finkelstein coordinates the 
proper coordinate transformation is given by 


y = t r log(r — 1) 
u = r. 


(25a) 

(25b) 


to consider additional coordinate systems. The angles d 
and if remain unchanged to keep the coordinates adapted 
to spherical symmetry. In the wave equation (21) partial 
derivatives with respect to r are meant that t is kept 
fixed, and vice versa. In the new coordinates, we have to 
work with partial derivatives with respect to u where v 
is kept fixed, and vice versa. Thus, we need the following 
transformation formulas for partial derivatives: 


d 


dv d 

du d 

Ft 

r=const. 

dt dv 

u=const. du 


/(r, i) I;+!)(■•. f)|;. 


(23a) 


The new time coordinate v is constant on ingoing radial 
null geodesics. With the transformation above the wave 
equation (|24|) reduces to the simpler form 


r — 1 


,(r-l) 92 


9^ 

dj .2 j. Qrdv 

(r — 1) d 


+ - 


dr 


- y{r) 


4'(z;,r)=0. (26) 


d 


dv d 


du d 


transformation 


dr 

t=const. 

dr dv 

u=const. 

' dr du 

i;=const. 

V = t-\- 2y/r log 

u = r. 

Fr 



=:f(r.t)E+g 

f P ^ 

(23b) 

Vr + 1 


Painleve-Gullstrand coordinates 

To consider Painleve-Gullstrand coordinates we have 
to use instead the slightly more complicated coordinate 


(27a) 

(27b) 


If, for the sake of readability, we omit the arguments of 
the functions / and the wave equation reads 


r — 1 
r 

r — 1 
r 

r — 1 
r 

r — 1 
r 

r — 1 
r 


f - /' 


/2 .2 
9-9 


f 9 ' - 2/5 


9^'I'(m, v) 
dv"^ 

9^4'(m, v) 
du^ 

d'^'^{u,v) 


dudv 


r + - / 


- 1 ) 


-9-9 


94'(u, v) 
dv 

d^{u, v) 
du 


—V{r)'S{u,v) = 0, 


(24) 


Calculating all necessary first and second derivatives we 
obtain 


r- ly d^ (r- 1) d^ (r-1) 
r ) dr'^ r dv'^ drdv 


{r — 1) d {r — 1) d 
2r^/'^ dv dr 


V{r) 


vl/(u, r) = 0 


(28) 


as the corresponding wave equation in PG coordinates. 
Its qualitative features are similar to those of the wave 
equation in EF coordinates. However, as the latter looks 
simpler we will concentrate in the following on the EF 
coordinates. 


Kruskal-Szekeres coordinates 

The transformation to Kruskal-Szekeres coordinates 
(in the outer region r > 1) is achieved by 


where r has to be considered as a function of the 
new coordinates, i.e., r = r{u,v). Now, we can in¬ 
troduce Eddington-Finkelstein, Painleve-Gullstrand and 
Kruskal-Szekeres coordinates by specifying the two func¬ 
tions f{r,t) and g{r,t), respectively. For an overview 
of such regular coordinate systems for the Schwarzschild 
spacetime the reader is referred to, e.g., the work in |18j 
and references therein. 


V = \/r — \ e”/^ sinh(f/2), (29a) 

u = y/r — 1 e”/^ cosh(t/2) (29b) 


and, thus, equation (24) simplifies to 
^2 




(r- 1) 


V du^ dv"^ 


4e ”F(r) 


4'(u,z;)=0. (30) 
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Here, r has to be considered as a function of u and u, Electromagnetic perturbations 

implicitly given by the relation 

To describe the propagation of electromagnetic waves 
(r — 1) e*" = M — V , (31) we consider the source-free Maxwell equations 


and we recognize the wave operator in double null coor¬ 
dinates d^jdu^ — d^ldv^. In the interior region (r < 1) 
u and V have to be interchanged and this corresponds 
to a change of signs in the wave operator. While in the 
four other coordinate systems the wave equation can be 
reduced with a separation ansatz to a stationary equa¬ 
tion by splitting off a time factor, this is not true in KS 
coordinates. For this reason, we will not consider the KS 
coordinates in the rest of the paper. 


C. Linear wave equations for scalar, 
electromagnetic and gravitational perturbations 


= 0 , 
= 0 . 


(35a) 

(35b) 


Owing to the conformal invariance of (35a I and (35b) we 
can use the metric 

,, ,, , , -|- sin^ 'dd(p‘^) , , 

= -dt^ -f drl + ^ , (36) 

(r - 1) 

which is conformally equivalent to the Schwarzschild met¬ 
ric. Following Stephani [^, the source-free Maxwell equa¬ 
tions can be reduced to the Debye-equation, which is a 
single differential equation for a scalar function called De¬ 
bye potential. For the background spacetime described 
by (36) this Debye equation is 


For each (massless) perturbation with any of the spin 
values s = 0,1, 2 a wave equation of the form 0 can 
be derived. These wave equations are valid in the sense 
that the back reaction on the spacetime geometry is ne¬ 
glected and we do consider the fields as small perturba¬ 
tions on a fixed background spacetime. In the follow¬ 
ing, we sketch how to derive the respective wave equa¬ 
tion on Schwarzschild background for a chosen spin value 
and combine the results in the end to obtain the general 
Regge-Wheeler equation that contains the spin as a pa¬ 
rameter. 


Scalar perturbations 


(^- 1 ) 
^3 


/ „ 9 9^ 1 9^ \ 

(^COt ^ ^ ^ j 


02 a 2 - 


n = 0 , 


(37) 


which is given in [12] as well. To obtain the general so¬ 
lution of Maxwell’s equations we have to construct the 
four-potential of the electromagnetic field [S] 

Ap = {xAvp - v''up)dyli + , (38) 


where H and are two independent solutions of the 
Debye-equation and (v^) = (0,1,0, 0), (u^) = (1, 0, 0,0). 
Thereupon, we can obtain the components of the electro¬ 
magnetic field as usual by 


For (massless) scalar perturbations the wave equation 
is a result of the Klein-Gordon-equation 


□$ = 


d 

dx^ 




d ' 
dx’' 


d> = 0, 
(32) 


where g is the determinant of the Schwarzschild metric. 
(Note that we use the signature convention 
If we expand the field into scalar spherical harmonics 
l)m('d, (/j) according to [5] 


OO I T \ 

$ = $(t,r*,'d,(/j) = ^ ^ (33) 

l—O m— — l 

and insert this ansatz into the Klein-Gordon equation, we 
notice that the remaining function 5'/(t,r*) has to fulfill 
the wave equation ( |17| ) with a potential given by 

U..oM=^(i(i + l) + )) - (34) 


Fni/ — dpAi, di/Ap 


(39) 


Since the Debye potential H in (371 is a scalar function, 
we can expand it into scalar spherical harmonics accord¬ 
ing to 


OO I 

n(t,r*,d,(p) = EE . (40) 

1 — 1 m— — l 


Inserting this ansatz into the Debye equation, we notice 
that the function 4'/(t, r*) has to fulfill the wave equation 
0 with a potential given by 

KM=i(r) = ^^^Z(Z + 1). (41) 


Gravitational perturbations 

In a linear perturbation approach, as considered by 
Regge and Wheeler [1], we write the total metric as 
the sum 


5/41/ - 5/44/ T dpi/ 
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of a background metric and a small perturbation 
We neglect all higher order terms in and its deriva¬ 
tives. A rather lengthy calculation, involving the expan¬ 
sion of the perturbation in scalar, vectorial and tensorial 
spherical harmonics and restricting to odd parity pertur¬ 
bations yields the wave equation ( |17[ ) with the potential 

= + (42) 

as a result of the Einstein vacuum field equation. The 
(odd-parity) components of the perturbation can be 
calculated in terms of the solutions to this wave equa¬ 
tion. A similar treatment for the even parity perturba¬ 
tions leads to the Zerilli equation [20] which, however will 
not be considered in the present paper. 


The general Regge-Wheeler equation 


Combining the previous results for different spins we 
obtain the general time-dependent Regge-Wheeler equa¬ 
tion 


drl 


dt^ 


Vsiir) 


= 0 , 


(43) 


where the spin-dependent potential is now given by 


UiM = lV^(i(i + i) + 


(44) 


The differential equation (431 is a wave equation of the 
form 0 and can, thus, be transformed to the corre¬ 
sponding wave equations in Schwarzschild, Eddington- 
Finkelstein, Painleve-Gullstrand and Kruskal-Szekeres 
coordinates using the results of section (IIB). 


III. RADIAL EQUATION IN TORTOISE AND 
SCHWARZSCHILD COORDINATES 


A. Stationary Regge-Wheeler equation 


A lot of the phenomenology can be understood by in¬ 
vestigating the wave equation in the simplest form, i.e., 
written in the tortoise coordinate. Hence, we shall briefly 
recap the properties of this differential equation. Assum¬ 
ing a harmonic time dependence 




(45) 


in (43) leads to the stationary form of the Regge-Wheeler 
equation (RWE) that is 


dr 2 


Vsi{r) 


Rujsl (^* ) 0 ■ 


(46) 


This radial equation is an ordinary, second order, linear 
and homogeneous differential equation in the form of a 



FIG. 1. The Regge-Wheeler potential as a function of the 
tortoise coordinate r* for I — 2 and spin s = 0 (dotted), s = 1 
(dashed) and s = 2 (solid). 



FIG. 2. The Regge-Wheeler potential as a function of the 
tortoise coordinate r* for s = 0 and I — 0 (dotted), I = 1 
(dashed) and I = 2 (solid). 


Schrodinger equation with “energy” The shape of 
the involved potential is shown in Fig. ([^ for all three 
spin values and a partial wave index ^ = 2 as an ex¬ 
ample. The maxima occur near the lightsphere radius 
r = 3/2 (r* « 0.81). While for s = 1 the maximum is al¬ 
ways located exactly at that radius, for s = 0 it is found 
slightly below and for s = 2 slightly above, but in any 
case it coincides with the lightsphere radius in the limit 
1 —^00. The value of the potential increases with the 
partial wave number (angular momentum index). Figure 
(H) shows this /-dependence of the maximum for s = 0. 
The potential decreases exponentially near the horizon 
(r* —> —oo) and behaves like for r* —>■ -|-oo. Hence, 
the asymptotic solutions in these limits are 

„ . 4 / r*^-oo 

r, ^+oo' 

This is, together with the time-dependent factor 
a superposition of in- and outgoing waves, where terms 
~ e”^’'* describe waves that propagate to smaller radii 
and terms ^ 0 *“”'* describe waves propagating to larger 






















radii with increasing time. To consider a black hole with¬ 
out quantum effects we have to assume causal boundary 
conditions to ensure that no waves emerge from the hori¬ 
zon. Thus, the asymptotic solution that satisfies this 
causal boundary condition is the IN-mode [5] 


RZir*) 


A i 

y^uji ' 


r* —)• —00 
r* —>■ -1-00 

(47) 


Usually, the amplitude at the horizon is normalized such 
that = 1 and we can define the reflection and 

transmission coefficients [3] 


B. Radial equation in Schwarzschild coordinate 


As a result of section II B| we know that the analog 
to the time-dependent RWE (coordinates t and r») in 
Schwarzschild coordinates {t and r) is obtained by using 
equation (21) together with the Regge-Wheeler potential 
(441. We can reduce this wave equation, again, to a radial 
equation with the ansatz 


as was done before to derive the stationary RWE. This 
leads to the radial equation 



1 

2 

A out 

2 

II 

N 3 

II 

3 

Ain 

) ^ujI = \R\ujI — 

'^ujI 

Ain 

^iol 

■ (48) 


dr^ dr 


The complex conjugate of (471 gives the OUT-mode, 
Considering the Wronskian 
we can derive the relation 


'^ujI + 1 


= 1 


(49) 


-(r- 1) 


/(; + !) 



Rusi{r) = 0 , 


(52) 


which is discussed in the literature for example by Leaver 
[3 [TO] and by Eiziev liin]. Leaver uses the ansatz 


that corresponds to flux conservation. To consider scat¬ 
tering at a black hole we have to find a solution of the 
RWE (46) that satisfies the causal boundary condition, 
i.e., that has the asymptotic behavior (47). Erom this 
solution we could read off the amplitudes A“ and 
to obtain the reflection and transmission coefficients or 
the phase shifts Si 


^2^Sl 


A out 
Ain 


(50) 


which give rise to the (total) scattering and absorption 
cross sections of the black hole. 

To obtain quasi-normal modes (QNM) and calculate 
their frequencies one usually restricts to purely ingoing 
boundary conditions: No waves emerge from the hori¬ 
zon, but no waves are coming in from spatial infinity 
as well. Hence, an initial perturbation radiates into the 
horizon and into spatial infinity and should decay in the 
cause of time in order to have a black hole that is sta¬ 
ble against the considered perturbation. This leads, of 
course, to imaginary frequencies, where the real part de¬ 
scribes an ordinary oscillation and the imaginary part 
(with the proper sign) assures the temporal decay. The 
proper asymptotic solution is then given by 


R.si{r) = r^+\r-l)-^^y^siir) 


(53) 


to transform equation (521 into a generalized spheroidal 
wave equation for j/,^s;(r)7which is nothing but a special 
form of the confluent Heun equation (CHE), even though 
Leaver never stated this nor did he use the term “con¬ 
fluent Heun equation” in his work. Eiziev on the other 
hand makes the ansatz 


Ru^si (r) = (r - (r) (54) 

and gives the confluent Heun equation in the Maple-form 
for the remaining function Hi^si(r) and its solutions in 
terms of the standard confluent Heun function. It should 
be stated that in both cases the ansatz is not well defined 
at the horizon (r = I) due to the term (r — 1)=*=“. In 
the following, we will briefly recall these known solutions 
since we can use them for a comparison of our results 
later. With the ansatz of Eiziev the function Hi^si{r) has 
to satisfy a confluent Heun equation in the Maple-form 
([^ with the five parameters 

a = 2iuj , b = 2s , c= 2ioj , (55) 

d=2j^, e = s^ — l{l + \). (56) 


R. 


QNM 

)sl 


(r*) = 




j^QNM 

ljI 




oo 


-foo 


(51) 


Unfortunately, no exact solution of the RWE (46) is 


known due to the involved potential (44). Thus, to solve 


the scattering or quasi-normal mode problem the only 
possibilities seem to be either approximation methods 
(e.g., WKB [TO], Born [^, ...) or numerical solutions. 
But yet another method is shown in the next sections 
where we derive exact analytic solutions, given in the 
form of convergent series, after changing the coordinate 
system and thus the underlying wave equation. 


The two regular singularities are located at the origin 
r = 0 and at the horizon r = 1, while the irregular sin¬ 
gularity is found at spatial infinity r = oo. We can give 
the GRS of this equation using the general result ([^ to 
complement the work in [8] and [9] . This GRS then reads 

/II 2 

0 1 oo 

0 0 s + 1 

—2s —2iuj s -I- 1 -I- 2iuj 

0 

V —2ioj 


( 57 ) 

























9 


Thus, the two local Frobenius solutions at the horizon 
can easily be constructed manually with the help of the 
above GRS or given in terms of the standard confluent 


Heun function as follows: 

1) = (r- l)"‘=HeunC(-a,-c,6,-d,e + d, 1 - r) 
1) = HeunC(—a, +c, 6, —d, e + d,l — r). 


Now, we use again (54) and obtain the full solutions of the 


Schwarzschild radial perturbation equation (52) around 


the regular singularity at the horizon. These solutions 


are 


1) = ^^^HeunC(—2ja;, —2ia;, 2s, —2a;^, 2a;^ + + 1), 1 — r), (59a) 

Rliiir-, 1) = r"+ie*‘^(’'+^°s(^-i»HeunC(-2ia;, +2iuj, 2s, - 20 ;^, 2u;^ + s^ - l{l + 1), 1 - r). (59b) 


It is important to remark that both solutions are not 
well defined exactly at the horizon. While the mod¬ 
ulus approaches one, the phase is not defined: In the 
complex plane infinitely many turns on the unit cir¬ 
cle are performed when approaching the singularity at 
r = 1. Thus, in these coordinates it is not possible to 
construct wavelike solutions that propagate through the 
horizon. Furthermore, the representation of the solutions 
in terms of the HeunC-function is only useful in the do¬ 
main I r — 11 < 1. 

The asymptotic Thome solutions at the irregular sin¬ 
gularity r = oo can be constructed, as shown in Sec. [T^ 
using the characteristic exponents of the second kind in 
the third column of the GRS (ISTl) 


HLiir-,oo) =r '^Pkr , 

k^O 

oo 

HiUr; oo) = ^ . 

k=0 


Thereupon, the full Thome solutions of the Schwarzschild 
radial equation at z = oo are 


OO 


RLiir;oo) 

Riiiir;oo) 


^ ^ ^ 

^lujr 

— Q-i‘^ir+log{r-l)-2log(r)) 





(6Ia) 

k^O 


oo 

(61b) 


which coincide with those given in |S]. These series are 
convergent in the limit r —>■ oo. For the expansion co¬ 
efficients pk and (Tfc recurrence relations can be derived 
and the normalizations po and ctoj be., the amplitudes at 
spatial infinity remain free to choose. 

Combined with the time-dependent factor the 

solutions and i?f,^;(r;oo) describe waves that 

propagate into the respective singularity, i.e., waves going 
into the horizon and waves going into spatial infinity. The 
solutions R^j^iir;!) and R^^;(r;oo) describe waves that 


-RLi{r;l) RLiir-,oo) -» 

RUsI (r; 1) Riil (r; oo) * - 



FIG. 3. In- and outgoing properties of the local solutions 
to the Schwarzschild and Eddington-Finkelstein radial equa¬ 
tions. 


emerge from the respective singular point. Figure (§ 
sketches these in- and outgoing properties as a mnemonic. 


IV. RADIAL EQUATION IN 
EDDINGTON-FINKELSTEIN COORDINATES 


To obtain the radial perturbation equation in 
Eddington-Finkelstein coordinates we use the wave equa¬ 
tion (26) together with the Regge-Wheeler potential (441. 
Thereupon, we perform the ansatz 


4'si(r;, r)=e R^siir) 


(62) 


to separate the Eddington-Finkelstein time coordinate v. 
Thus, the remaining function Rujsiif) has to fulfill the 
differential equation 


+ 

/'I — 2iojr^ ^ 


dj.2 

V r{r - 1) y 

dr 


l{l -1-1) -I- (1 - s^)/r 
r{r - I) 


Rijsi{r) = 0 . 


(63) 


This radial equation exhibits two regular singularities 
which are located at r = 0 and r = 1, since these points 
are poles of first order in the second term and poles of 
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at most second order in the last term. By means of the 
substitution r i—> 1/Q we obtain 


'd^ f‘iC,^ - 2C, - 2iuj\ d 

C"(C- 1 ) ) dC 

/ ;(/ + l) + C(l-s^) Y 

A C^(C-i) A 


RiosiiC) 


0 (64) 


which yields a pole of second order at C = 0 in the co¬ 
efficient of the first derivative term. Hence, we conclude 
that the original differential equation (63) possesses an 
irregular singularity at r = oo. With two regular singu¬ 
larities at r = 0,1 and the irregular singularity at r = oo 
the radial equation in Eddington-Finkelstein coordinates 
belongs to the class of singly confluent Heun equations 
as well. It is indeed possible to find a substitution that 


transforms (63) into a standard form of the CHE. One 
such possibility simply is 


Riusi{r) = ( 65 ) 

and yields a CHE in the form ([^ with the set of param¬ 
eter 


a = (s -I- 1), /3 = 2iuj, 7 = 1 -I- 2s , 

(5=1 — 2iu }, q = s(s -|- 1) — l(l + 1). (66) 


A. Local solutions 


Using the five parameters in (66) we can immediately 


give the corresponding GRS for the differential equation 


following the general result 

/II 2 \ 

01 oo ;r 

0 0 s-f 1 

—2s 2iuj s -I- 1 — 2iuj 

0 

V 2iijj / 


(67) 


This GRS already reveals a lot of useful information: 
Since s can only take one of the integer values 0,1, 2 the 
two indicial exponents at the regular singularity r = 0 
differ only by an integer. Hence, both local Erobenius 
solutions constructed at this point are not linearly inde¬ 
pendent anymore |16j . Moreover, we can already see that 
it will be possible to construct a solution that is regular 
at the horizon. This is due to the indicial exponent 0 
in the second column of the GRS and the transforma¬ 
tion (65) that is regular at the horizon r = 1, unlike in 
the Schwarzschild coordinate case before. It is also pos¬ 
sible to give the radial equation in EE coordinates in the 
Maple form of the GHE. In this case the parameter for 
the confluent Heun equation in the form (§ become 


a = —2iu}, b = 2s, c = —2iuj, 
d = 2a;^, e = s^ — l{l + 1). (68) 

Using these five parameters we can give the solutions 
around the regular singularity at the horizon in terms of 
the standard confluent Heun function as follows: 


Risiir-, 1) = yisiir-, 1) = HeunG(25a;, -2iuj, 2s, -2a;^ - l{l -k 1) -f 2w^ 1 - r), (69a) 

Riiiir; 1) = yiiiir; 1) = r^+\r - 1)^“ HeunG(2*a;, 2*cc, 2s, -2^^, s^ - + 1) + 2lo^, 1 - t) . (69b) 


Both solutions are independent and differ remarkably in 
their behavior at the horizon: While the first solution 
is well behaved and takes simply a constant value, the 
second solution A^;(r; 1) performs infinitely many turns 
on the unit circle in the complex plane and has no well 
defined phase when approaching r = 1. Figure ([^ shows 
the real part and the absolute value of both solutions for 
some arbitrarily chosen values of the parameter w, s, 1. 


Around the irregular singularity at spatial infinity we 
can construct the Thome solutions, using the character¬ 
istic exponents of second kind in the third column of the 


GRS (67). The result is 


RLiir-, oo) = £ p,(a;, s, 1) r-^ , (70a) 




Rusiir-,oo) = 


'^ak{uj,s,l) 




(70b) 


fc =0 


and the coefficients pk and au can be calculated using 
recurrence relations obtained from the differential equa¬ 
tion, again. 

We labeled all the solutions such that we get the same 
scheme ([^ as in the Schwarzschild case before: The reg¬ 
ular solution at the horizon As;(r;l) describes waves 
that propagate into the horizon, while the other solu¬ 
tion RUslir-, 1) is associated with waves that emerge from 
the horizon and propagate to larger radii. At the ir- 















11 




r 


FIG. 4. The two local solutions around the horizon for 
uj = 30 ,1 = 30, s = 0: the modulus (solid) and the real part 
(dotted) are shown. 


regular singularity at spatial infinity the two solutions 
{r; oo) describe in- and outgoing waves as well: The 
first solution i?j^^;(r;oo) describes waves that propagate 
into the singularity (outgoing) and the second solution 
RU^i{r;oo) is associated with waves that emerge from 
spatial infinity and propagate to smaller radii (ingoing). 


and Slavyanov and Lay m gave a general guide for ap¬ 
plying the method of Jaffe to any kind of Heun equation, 
including all the confluent cases, and they worked out 
what is called the central two point connection problem 
(CTCP) in great detail. The derivations in this section 
are based on the work of Slavyanov and Lay and we will 
put our radial equation in such a form that the results in 
usi can be used. 

The question of interest is now: is it possible to “con¬ 
nect” the local solutions at the horizon and those at spa¬ 
tial infinity and obtain expressions that are valid on a 
larger domain or cover even the entire positive real axis 
as the domain of convergence? The basic idea is to “glue” 
Frobenius- and Thome-type solutions together. In this 
way we encounter a special case of the central two point 
connection problem. To solve this problem we will apply 
the following five step-procedure: 

1. Shift the regular singularity and reorder the differ¬ 
ential equation. 

2. Perform an s-homotopic transformation of the de¬ 
pendent variable that involves the indicial exponent 
of the respective Frobenius solution at the horizon 
and the characteristic Thome exponents (of order 
zero and one) at spatial infinity. 

3. Mobius-transform the independent variable such 
that the regular singularity remains unchanged at 
the origin, but the irregular singularity is moved to 
the boundary of the unit circle. 

4. Solve the resulting differential equation using a se¬ 
ries expansion around the origin that converges on 
the open unit circle. 

5. Construct, thereupon, the total solution of the orig¬ 
inal radial equation by “inverting” the previous 
transformations. 


B. Analytic continuation of the local solutions 


It is quite unp leasa nt th at the two local solutions 
RI^Ii {r\ 1) in (69a) and (69b) have such a small radius of 
convergence which makes them useful only in the region 
|1 —r| < 1. This is a consequence of the fact that the only 
known representation of the confluent Heun function is 
the one in terms of a Frobenius series. When we consider 
scattering of radiation at a black hole or when we are 
interested in its quasi-normal modes we have to impose 
boundary conditions at spatial infinity, or at least at a 
sufficiently large radius, and therefore the representation 
of the confluent Heun function HeunC as a Frobenius 
series is not appropriate. 

In this section we will work out a procedure that is 
related to a similar problem considered by Jaffe [22]. He 
calculated the bound states and energy eigenvalues of the 
hydrogen molecule ion. Actually, the same Idea was ap¬ 
plied by Leaver [ziiin] to obtain a continued fraction ex¬ 
pression for quasi-normal modes of black holes. Lay |23j 


1. Shift the singularities and reorder the differential 
equation 


We use the transformation r i—)■ z = r — 1 that maps 
the regular singularity from r = 1 to z = 0, while the 
irregular singularity remains at z = oo. Here it is impor¬ 
tant that no other singularity is located between those 
two on the positive real axis. This yields the differential 
equation 


dz^ 




■ z* 



_d 

d^ 



yuisi{z) = 0 


with the parameters 


(71) 


— 1, Gq — — ‘liuj , — 1 — ‘lioj , ^2 — 1-4“ 2s , 

Cl = s(s + 1) - l{l + 1) - 2iuj{s + 1), (72) 

C 2 = 1(1 H“ 1) — s(s + 1 ), 
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Frobenius solutions 


KJ’-A) /'i= 0 -; 7 - 


fi2 = “(s + 1) + 2iuj 


Thome solutions 
-p.v = 2iu} i?i,,(r;oo) 








1 ) = 2 ia; . 


fi 2 = —(s + 1) — 2iui 


-^v = {) R".i(r-,ao) 


FIG. 5. A scheme showing the proper choice of the expo¬ 
nents V, fii and p2 related to local solutions for the connection 
problem. 


and is done to match the conventions used in [15] . 


new differential equation for ^{x), that is 


x^ix-l)‘ 




lix) 


dx^ 


^ ^ ^kX 


k=l 


k d- ^ojsiix) 

dx 


+ J2^kX^^u.si{x)=0, (75) 


k=l 


where the coefficients flk and are given by 


fli = 1 — 2iu! + 2/ri, 

112 = 2(s -I- i/ -I- /i 2 - Ail - 1), 

Ha = 1 — 2 s — 2 fi 2 , 

Aa = s(s + 1 ) — l{l -l- 1 ) 2 (s -|- Ai 2 )(/ti — k') 

- -I- Ail + 1^2 

A 2 = At2(2s-I-Ai 2 ) ■ (76) 


2. S-homotopic transformation of the dependent variable 


4- A series expansion ansatz 


We perform an s-homotopic transformation of the de¬ 
pendent variable yu>si{z) that contains indicial exponents 
at 2 ; = 0 and characteristic exponents of the second kind 
(order zero and order one) at 2 ; = 00 . This transforma¬ 
tion is given by 

yujsl{z) Eujsl{z) 

yujsi{z) = e‘'^ z'"! (z -f 1)'"" S^s/(z), (73) 


We use a series expansion ansatz for the unknown func¬ 
tion S^s/(x) 

00 

E^siix) ='^^k{oJ,s,l) x'' (77) 

fc =0 

and insert this ansatz into the differential equation. 
Thus, we obtain a three term recurrence relation for the 
series coefficients ^k 


where the exponents v, /ii, y .2 have the following meaning: 


Hi is the indicial exponent of the chosen Frobenius 
solution at z = 0. Thus, we can have either fxi = 0 
or Ail = 2ioj according to the GRS (67). 


- ir is the characteristic Thome exponent of order one 
and depends on the chosen Thome solution at z = 
00 . Hence, we can have either 1 / = 0 or 1 / = 2iu!. 


- For large z we get z^i(z -I- 1)^= « 
have to adjust the sum Ati +Ai 2 to be the character¬ 
istic exponent of order zero of the chosen Thome 
solution. Hence, the possible values for fj ,2 are 

Ai 2 = — (s + 1 ) or At 2 = — (s + 1 ) ± 2 iw . 


The scheme in Fig. shows which combination of the 
three exponents must be taken for the desired choice of 
Frobenius and Thome solutions to be connected. 


3. Mobius transformation of the independent variable 


Next, we apply the Mobius transformation 


z X 


z 

z + 1 


(74) 


that maps the irregular singularity from z = 00 to x = 1 
and leaves the regular singularity at x = 0. This yields a 


^0 = arbitrary 
OoCi + Po^o = 0 

Ofc^fe-l-l + Pk^k + Ikik-i = 0 , (78) 

where \/k > 0 we have 


ao — Hi -1-1, /3o — Ai, 


CTfc — 1 + 


Pk — ~2 -|- 


Hi -f 1 

k 

H 2 -l- 2 


Hi 

F’ 

+ — 
fc 2 ’ 


7fc = 1 - 


H 3 — 3 A 2 — H 3 -|- 2 


(79) 


The solution ([T^ with the recurrence relation ( |78| ) is ab¬ 
solutely convergent on the open unit circle around x = 0 , 
i.e., the radius of convergence is the distance between the 
two considered singularities at x = 0 and x = 1. Choose 
any x with |x| < 1 , then 


lim 

k—>-OQ 


^kX^ 


= |x| < 1 


and convergence for |x| < 1 is assured |7]. However, uni¬ 
form convergence will strongly depend on the properties 
of the coefficients ^^(w, s, Z), since 


lim 

lim 

C/c+lx''+^ 

= lim 

^fc+i 

k—xx) 

fC —>-1 

ikx^ 

k—^oo 
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Hence, the solution is uniformly convergent if 
s,l) exists and is finite, which will depend on 
the values of w, s and 1. The sum will be finite and, 
thus, uniform convergence will be assured if s, 1) is 
a minimal solution of the recurrence relation (78) [see, 

e.g.,mm- 


5. Final solution of the radial equation 


Finally, we can construct the total solution to the ra¬ 
dial equation in EF coordinates (63). We have to use 


V. APPLICATIONS 
A. Black hole scattering 


As an example, we consider the scattering of a mass¬ 
less scalar field at a Schwarzschild black hole. Numerical 
studies of such scattering processes can be found for ex¬ 
ample in [23] and |2S| . In analogy to the usual scattering 
theory used in quantum mechanics and nuclear physics 
we consider an ingoing plane wave coming from spatial 
infinity as initial condition. Such a plane scalar wave is 
in Schwarzschild coordinates for r —>■ cx) given by |26] 


r — 1 


X = 


1 


(80) 


to recover the original radius variable r. Thereupon, we 
obtain 




(r - 1)'^ 




(81) 


The exponents v, , /i2 must be chosen according to the 
scheme in Fig. and the requirements of the physical 
situation that is to be described. The solution (81) is 


absolutely convergent on the open interval r g]0.5, (X)[ 
owing to the absolute convergence of E!(xHor |a;| < 1 and 
the mapping x —>■ r a.s explained in Fig. (|m. In the entire 


convergence radius of Frobenius solutions {r] 1) 



|r- 1| < 1 

1 


1 0 

0.5 r = l 

1 

■2 

CX) 

1 — X 

\ \ 



r — 1 

\ 



r 

-li 

1 

z = 0 

( 

1 


1 

1 

|z|<l 

1 

1 


convergence radius of E.^si{x) 


FIG. 6. A scheme showing the mapping from a; to r and the 
solutions Ruisiir) and Suisiix) together with their domains of 
convergence. 


region where the Frobenius solutions and the analytically 
continued solutions overlap, they do coincide. Hence, we 
can interpret our solutions (81) as new representations of 
the confluent Heun function HeunC with larger domain 
of convergence. In the next section we show how to use 
these solutions for different physical applications. 


$ 


plane 


T 

1=0 


l\2l + l) . f l7T\ 

-sm cur^ —— Fiicosv ). 

Lur \ 2 J 

(82) 


This is a superposition of partial waves with angular mo¬ 
mentum index 1. Each partial wave will be reflected 
and/or absorbed by the black hole in a different manner 
and this will change the composition of the total wave 
field since different modes will be modified with different 
amplitudes that refer to scattering and absorption coef¬ 
ficients. This can be best understood by regarding the 
RWE (46) and the involved potential. For a given fre¬ 


quency u} of the incoming wave the potential depends on 
the wave index I as shown in Fig. (H). As long as the 
energy of the incoming wave is larger than the poten¬ 
tial peak for a given Z-mode, this mode will be (mostly) 
absorbed by the black hole. When the energy becomes 
smaller than the potential peak for a fixed Z-mode, this 
mode will be (mostly) reflected by the potential barrier. 
For the special case we expect the transmis¬ 
sion and reflection probabilities to be « 0.5. 

To derive the transmission and absorption probabilities 
we must use the solution (81) with the exponents chosen 
to be 


= 1 ^ 1 = 0 , fi2 = -is+l) 


/ 1 \ 

i?u.si(r) = ^^fe(w,s,Z)l^^—^— (83) 


fc =0 


to match the causal boundary condition, i.e., to ensure 
purely ingoing waves at the horizon. The full wave field 
in EF coordinates is then given by 

$^e-“’^f;^^^i?„«;(r)Pi(cos^5). (84) 

i=0 

We can, of course, transform this solution to usual 
Schwarzschild coordinates simply by using v = t + r^f = 
t + r + log(r — 1). This might be done since the plane 
wave is given in Schwarzschild coordinates as well. (The 
other option is to give the plane wave in EF coordinates.) 
















14 


Thus, we have 

1=0 ^ 

(85) 

and we know that in the limit r —> oo the asymptotic 
form must be 

$ ^ ^ (2[^ P;(cos^»). 

/=o ^ 

( 86 ) 

Now, we can proceed as in usual quantum mechanical 
scattering theory: At large distances the total wave field 
shall be the sum of the incoming plane wave and a back- 
scattered part due to reflection at the black hole 

giwr, 

^ ^piane + e"*‘^‘/(^?)- for T ^ OO , (87) 

r 

and we define the complex valued phase shifts Si by 

^ - <I>pia„e ~ e-“*e“’-* f2 - l) Pz(cos^?). 

( 88 ) 


We can show that this implies 


2i5i 


= (-l) 


l-\-l 


A out 

Ain 

^ujI 


(89) 


Hence, we have to calculate the transmission and reflec¬ 
tion amplitudes to obtain the phase shifts. We 

can do so by matching our solution (851 with the radial 


function given by (831 to the asymptotic form ( 861 . Even 


though the radial function in the form (831 is not in gen¬ 
eral convergent in the limit r —> oo for an arbitrary choice 
of the parameter to, s and I we can fit it to the asymptotic 
behavior on an interval for large values of r and for fixed 
mode number I and frequency w. With this procedure we 
can calculate the phase shifts in a semi-analytical way 
for all desired combinations of the parameters. Using 
the phase shifts or the reflection and transmission am¬ 
plitudes we can calculate the transmission and reflection 
coefficients 


T 


col — 


\T\ 


2 _ 
LOI 




2 


^col — 


\R\ 


2 _ 
Cdl 


A out 


2 


A 


in 

cul 


They describe the probability of being transmitted 
through the potential barrier, i.e., being absorbed by the 
black hole and the probability of being reflected back 
to spatial infinity. The transmission coefficients are 
usually called graybody factors, since they are important 
when we derive the deviations from a thermal spectrum 
of Hawking radiation. Figure 0 shows the calculated 
graybody factors as functions of the frequency for differ¬ 
ent wave modes. Using the graybody factors we can, for 



FIG. 7. Graybody factors as functions of the frequency for 
/ = 0 (thick solid) up to 1 = 5 (thin dashed). 


example, calculate the total absorption cross section as a 
sum over all partial contributions 

+ (90) 

In Fig. (|^ we show the partial absorption cross sections 
for different values of the wave mode starting from I = 0 
(on the left) to Z = 5. The results are normalized by the 
horizon area. The figure further shows the total absorp¬ 
tion cross section as a sum over all partial contributions. 
The result is normalized by the geometrical optics value 
of the cross section 277rM^, which is recovered in the 
limit of high frequencies as it should be. For small fre¬ 
quencies ojM —)■ 0 the partial absorption cross section for 
the I = s = 0 mode (lowest mode) approaches the horizon 
area of the black hole and therefore overcomes the prob¬ 
lem shown in Figure (7) in the phase-integral study by 
Andersson [26]. For higher frequencies our results match 
those given by Andersson and our results in Fig. ® are 
in agreement with those shown in Figures (2) and (3) in 
the early work by Sanchez 


B. Quasi-normal modes 


We can use the radial solutions (811 to obtain the 


quasi-normal mode frequencies of a Schwarzschild black 
hole for a spin s perturbation. To ensure the proper 
QNM boundary conditions we have to choose, according 
to the scheme in Fig. ([^, the exponents to be 


V = 2iuj, /ii = 0, /r 2 = — (s -I- 1) -I- 2iuj 


R' 


(r) = 


^Cfc(w,s, Z) ^ 
k=0 


(91) 


and we further demand the solution to converge in the 
limit r —>■ oo. The convergence property is fulfilled if 
the series coefficients are minimal solutions of the re¬ 
currence relation (78). This immediately yields the con¬ 


tinued fraction equation that was given by Leaver m 
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FIG. 8. Top: partial absorption cross sections normalized 
by the horizon area Ah for wave modes starting from I = 0 
(left) up to 1 = 5; Bottom: total absorption cross section as 
a sum over all contributions normalized by the geometrical 
optics value 27nM^. 


and the QNM frequencies ojqnm can be calculated as 
complex roots of this equation. Hence, we see that the 
method of Leaver fits perfectly into our framework of 
continued local Frobenius solutions. For details on quasi¬ 
normal mode frequencies the reader is referred to Refs. 
PSI [53] and [301 ^"^cl to the living reviews (OB [32] . 


C. Hawking radiation 


Using the method of Damour and Ruffini m we can 
calculate the spectrum of the Hawking radiation emitted 
by a Schwarzschild black hole. Now, we have to consider 
that part of the solution to the radial equation which was 
always neglected up to now: We need to use the second 
solution {r; 1) that is not regular at the horizon. This 
solution is given by (69b) and describes purely outgoing 


waves at the horizon r = 1: This is the outgoing Hawk¬ 
ing radiation that is obtained here as the solution of a 
classical wave equation on curved background geometry. 
Close to the horizon this solution behaves like 


- {r-iy 


(92) 


Following Damour and Ruffini, there exists a unique con¬ 
tinuation of this solution for radii r < 1, which yields 


r > 1 




(93) 


This continuation then describes an antiparticle of nega¬ 
tive energy inside going down the black hole and a par¬ 
ticle outside emerging from the horizon. The probabil¬ 
ity for transmission through the horizon is given by the 
squared modulus of the relative amplitude 

(= . (94) 

Using Feynman’s method [33] of scattering amplitudes 
one can deduce that this is also the probability for 
antiparticle-particle pair production just outside r = 1. 
In |34j Sannan showed how to get the mean number of 
particles that are radiated to spatial infinity in a given 
mode from the pair production probability p^j- The re¬ 
sult is 


N^i = 


qSttujM ^ ^ ’ 


(95) 


where the -I- is for fermionic particles and the — must be 
used for bosons. The bosonic result is exactly Hawking’s 
result [35] that was obtained using calculations involving 
the gravitational collapse. When we consider bosons, the 
spectrum is thermal except for the modihcation caused 
by the transmission probability (or graybody factor) , 
which causes an amplitude Hltering of different modes. 
Hence, the total spectrum is not really thermal anymore 
but is usually called graybody radiation (instead of black- 
body radiation with thermal spectrum). From ( [95| ) we 
can derive a power spectrum of the emitted radiation. 
Each massless boson takes away an energy of uj from the 
black hole and the change of energy per unit of time and 
frequency interval is minus the radiated power which is 
given by 


dE 

dujdt 


1 

27r 




(96) 


Here, p denotes additional degrees of freedom like the po¬ 
larization of the particle |33| ■ It should be mentioned that 
this radiation spectrum is valid for a time in which we 
consider the background geometry to remain unchanged, 
i.e., the mass of the black hole is taken to be constant 
during this time interval. Then, we may iterate and con¬ 
sider the mass to be decreased due to the evaporated 
energy and repeat the calculation of the radiated power 
with the decreased mass. In this way one can derive the 
mass of the black hole as a function of time and calculate 
the life time of black holes with given initial mass. Figure 
([^ shows the obtained power spectrum according to (96) 
for massless bosons with spin s G {0,1,2}. The results 
in Fig. (|^ match previous results, for example those ob¬ 
tained by Page [36] using purely numerical calculations. 
The major contributions in the spectra in Fig. © are 
due to the lowest I = s modes for each spin, respectively. 
The little “bumps” that we can see at higher frequencies 
are due to higher mode contributions. For a fixed spin 
the maxima of higher mode contributions are shifted to 
higher frequencies but their peak values decrease rapidly 
with increasing mode number. 
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FIG. 9. Power spectrum of Hawking radiation emitted by a 
Schwarzschild black hole for massless bosons with spin s = 0 
(dashed), s = 1 (dotted), s = 2 (dash-dotted) compared to a 
thermal spectrum for a cross section of 27'kM^ (solid). 


somewhat more involved. These solutions are given in 
terms of the standard confluent Heun function around 
the black hole horizon which is defined in terms of 
Frobenius series. As these series converge on the unit 
sphere in the complex plane centered at r = 2M, they 
cannot be evaluated at larger radii. We have shown how 
this problem can be overcome by using an analytical 
continuation procedure. This gave us analytical solu¬ 
tions in terms of convergent series on the entire domain 
up to (but not in general including) infinity. Based on 
these exact solutions we have re-obtained, in a semi- 
analytic fashion, several results on black-hole scattering, 
quasi-normal modes and Hawking radiation that had 
been found numerically earlier. In a follow-up paper 
we consider a massive scalar field in EF coordinates 
which allows us to define purely ingoing solutions that 
penetrate the horizon. Such a massive scalar field is a 
possible dark matter candidate and analytical solutions 
might be of interest in this area. In the same follow-up 
paper we further re-consider massless fields of spin 
s = 0, 1 or 2 but switch to the PG coordinates which 
are associated with infalling observers and therefore, at 
least in principle, with possible measurements. 


CONCLUSIONS 


The Regge-Wheeler equation is usually written in 
terms of the tortoise coordinate. However, it can 
also be written in other coordinate systems, such as 
usual Schwarzschild, Eddington-Finkelstein, Painleve- 
Gullstrand or Kruskal-Szekeres coordinates. While it is 
not possible to obtain regular solutions at the black-hole 
horizon in Schwarzschild coordinates (see the comment 
following Eq.(54)), we have shown that such solutions 
can be derived in EF coordinates. This is also possible 
in PG coordinates where, however, the expressions are 
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